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We show that given one O(a) improvement constant, 6 m , all the remaining quantities needed to define the 
renormalized and 0(a) improved dimension-3 quark bilinears can be obtained by studying the matrix elements 
of these operators between external quark states in a fixed gauge. 



1. INTRODUCTION 

One of the leading uncertainties in lattice cal- 
culations involves the connection between the lat- 
tice and the continuum renormalized operators. 
Current estimates jjj show that one-loop pertur- 
bation theory for 0(a) improved Wilson fermions 
underestimates quantities like Z P /Zg by ~ 10% 
at lattice scales of 2-4 GeV. Furthermore, pertur- 
bative estimates of the 0(a) improvement coeffi- 
cients are significantly different from their non- 
perturbative estimates. 

All the scale independent renormalization con- 
stants and the improvement coefficients in the 
quenched theory have been determined by impos- 
ing vector and axial ward identities on the lattice. 
This Ward identity method is computationally in- 
tensive and alternate methods are therefore desir- 
able, especially to determine the scale dependent 
renormalization constants (Zj, and Z P or Zg). 

An alternate well known method that gives all 
the renormalization constants Jpj and some of 
the improvement constants HfiJ involves calcu- 
lating the matrix elements of the quark bilinears 
between external quark states in a fixed gauge. 
This method is far more tractable computation- 
ally and the generalization to 4-fermion opera- 
tors is straightforward. Here, we show how this 
method gives all but one (b m as discussed in sec- 
tion 3.1) of the 0{a) renormalization constants. 
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2. METHOD 

We start by defining the notation. 0(a) im- 
provement of the theory is achieved by improving 
the action and the operators simultaneously. The 
improved renormalized quark fields, ip can be re- 
lated to the lattice quark field ip by 

i> = ^ 1/2 (i-^ am i) x 

[1 - ac^ip+mi) ~ ac NGI </)}ip , (1) 

where mj is an 0(a) improved quark mass. Thus, 
apart from a mass dependent renormalization 
constant, one needs (i) an equation of motion cor- 
rection, cL, that does not affect position space 
correlation functions at finite separations, and 
(ii) mixing with a gauge non-invariant operator, 
cngii that appears because the calculation is per- 
formed in a fixed gauge, [gj 

We also write the renormalized propagator as 

(T[W>\) = ;5 j - = -i*U* + *2 , (2) 

where T is the time-ordering symbol, and is 
the four momentum of the quark. Neglecting log- 
arithmic and higher order corrections in p~ 2 , the 
improved renormalized propagator at high mo- 
menta is given by 
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where Z m is the renormalization constant of m/. 
The important point to note in these expressions 
is that chiral symmetry prevents terms propor- 
tional to mj and m/ 2 in Ei and £2 respectively. 
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Figure 1. Plot of 02 versus p 2 before (diamonds) Figure 2. Plot of a 2 versus p~ 2 after subtraction 

and after (squares) subtraction of 0(a 2 p 2 ) arte- of 0(a 2 p 2 ) artefact for n — 0.1344. 

fact for k = 0.1344. 



Expanding the lattice quark fields in terms of 
the continuum field, lattice propagator is 

(T[V>?]) = -iaij>+a 2 (5) 
= Zjf,(l + 2b t f,am)(l + 2iacNGii>) x 
(-iaitf + (T 2 ) + 2ac^ . (6) 



From this the unknown constants Z, 



cngIi an d can be extracted as follows. We 
first expand <j\ and a 2 at large p 2 as 
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where the terms dropped are yet higher order in 



p 2 . From Eqns. [|-[|, we note that these lead 



ing and next to leading coefficients, a[ 
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the expansion of ci,2 in p 2 have the following 
dependence on the quark mass: 
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where, omitting terms of 0(a 2 ) 
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— 2aZ^,(b^ — cngi) 

= -2Z4f3 1 -2ac NG i(3' 2 ) 
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= -VaZ^b* - c Naj (20x - /%)](16) 
= -Z v ,(l + ai ) (17) 
= -aZ^ [6^ (1 + c NG i) - 

2 CA r G7 (l + 2ai -a 2 )] (18) 
2 0,0 = 2aZ m {Z^c N ai + c'^) (19) 
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2aZ m Z^,[b 1 p + cjvg/(1 + £*i)] • (22) 



The five renormalization and improvement con- 
stants can now be obtained once a\ is eliminated 

LOft LO.l NLOft LOft NLOA ■, 

using a l ■ , <t 1 , <r 1 , cr 9 , a, , and 



7 u 2 ? 2 



cr^ 1 ' 2 • Once the improved propagator is known, 
the remaining constants are obtained as discussed 
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3. EXAMPLE IMPLEMENTATION 

3.1. Lattice details 

To illustrate the feasibility of this method, we 
present a preliminary analysis of 50 32 4 quenched 
configurations at (3 — 6.2. We take csw — 1-614 
from the ALPHA collaboration j|] to define the 
0(a) improved action. Simulations are done at 
seven values of k = 0.131, 0.1321, 0.1333, 0.1339, 
0.1344, 0.1348, and 0.1350. The critical value of 
the hopping constant, k c = 0.135899, and the 
fourth root of the plaquette, uo = 0.88510 are 
taken from pj]. For the O(a) improved definition 
of quark mass we take 



mj = In 
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since the value of b m is close to its tadpole im- 
proved tree level value |J. We scale all lattice 
fermion fields by ^/2huq. The lattice momentum 
components are p = sin(27r//32) and we average 
over momentum combinations equivalent under 
the hypercubic lattice symmetry group. Our fits 
use the seventy distinct momentum combinations 
with j < 4. 

3.2. 0(a 2 ) subtraction and determination 
of the constants 

According to Eqs. [7] and ||, p 2 o\ and a-i are 
supposed to go to constants at large p 2 , however, 
the data, exemplified in Figure [l], show a linear 
behavior in a 2 p 2 at large p 2 . These 0(a 2 ) terms 
are removed by fitting to the large momentum 
behavior. 

To extract the desired coefficients cr^ 1 " we 
now fit these subtracted p 2 a\ and 02 as a function 
of p~ 2 (see Figure 2 for an example). The data 
show that one needs to keep at least terms up to 
p~ 4 in Eqs. f|| to obtain a reasonable fit. The 
resulting intercepts and slopes obtained are then 
fit against mi (e.g., see Figure 3) to obtain the 
various expressions defined in Eqns. [l5|-p2"[ From 
these fits, we find 
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Figure 3. Slope of 0(a 2 )-corrected 01 with re- 
spect to p~ 2 plotted versus mi. 



^ = 0.27(4) 
DISCUSSION 
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We have shown that once b m , or equivalently, 
mi is known, all other 0(a) improvement con- 
stants needed to define the rcnormalized propaga- 
tor can be determined. In previous work S, the 
constant cngi was left undetermined. Using per- 
turbation theory Sharpe Q has shown that the 
effect of this term is small. We show that this con- 
stant can be determined non-perturbatively and 
its value is indeed small. 
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